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AN ANALISIS OF THE EFFECTS OF AEEOEIASTICITY ON STATIC, 
LONGITUDINAL STABILITY AND CONTROL OF 
A SWEPT^ACK-WING AIRPIANE 
By Richard B. Skoog 

SUMMARY 


A theoretical analysis has heen made of the effects of aeroelastic— 
ity on the stick— fixed static longitudinal stability and elevator angle 
required for balance of an airplane. The analysis is based on the famil- 
iar stability equation expressing the contribution of wing and tail to 
longitvidinal stability. Effects of wing, tail, and fuselage flexibility 
are considered. Calctilated effects are shown for a swept-wing bomber of 
relatively high flexibility. 

Although large changes in stability due to certain parameters are 
indicated for the example airplane, the over-all stability change was 
quite small, compared to the individual effects, due to the counter- 
balancing of wing and tail contributions. The effect of flexibility on 
longitudinal control for the example airplane was found to be of little 
real importance in Ig flight, althoiagh in turning flight the effect was 
found to be commensurate with the stability loss. 


INTRODUCTION 

In the past, airplane configurations and operating speeds have been 
such that prediction of longitudinal— stability and —control characteris- 
tics usually could be handled without regard to aeroelastic effects. 

With hi^er flight speeds and the use of swept wings, aeroelastic effects 
are becoming sufficiently important to pose major problems in airplane 
design. Althou^ flexibility of swept wings has introduced the major 
problems to date, the flexibility of fuselage and swept tail surfaces 
may introduce problems approaching equal importance as speeds continue 
to increase. Most of the published work on predicting these effects 
appears to have been done by the British (references 1, 2, 3^ and 5)» 
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The object of the present study was to determine the mgnitude of 
the stability loss likely to be encoimtered on swept-wing aircraft of 
conventional configuration. To gain this end, it was decided to evaluate 
the aeroelastic effects for an airplane of relatively high flexibility. 

It was decided, further, to employ the simplest possible method of sta- 
bility analysis in order to gain a maximum physical appreciation of the 
factors involved in the net stability change for the airplane. 


The results of the afore-mentioned study are presented in this report 
together with the method of analysis en^iloyed. The net stability change 
is shown together with the individual contributions due to flexibility of 
wing, tail, and fuselage, both including and neglecting the effect of 
inertial loads. The method of analysis is based on the familiar stability 
equation expressing the contribution of wing and tail to longitudinal 
stability. The reader interested solely in the calculated results can 
turn directly to the section titled "Application to Example Airplane." 


A 


b 

c 

c 

av 




^t 


NOTATION 


wing aspect ratio 



ratio of net aerodynamic force along the airplane Z axis 
(positive when directed upward) to the weight of the air- 
plane 

wing span measured normal to plane of symmetry, feet 
section chord parallel to the plane of symmetry, feet 
average section— chord parallel to the plane of symmetry, feet 

mean aerodsmamic chord .of wing 
section lift coefficient 




1, feet 


c dy 


lift coefficient 


lift 

qS 


tail contribution to lift coefficient 


wing lift-curve slope, per radian 


t all lift 
‘It^t 
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.tail lift-curve slope, per radian 


pitching-moment coefficient, positive for nose-^ip moment 
^ moment A 
\ qSc / 

pitcMng-moment coefficient at zero lift 

rate of change of pitching-moment coefficient with lift 
coefficient 


H 




tail incidence [i.j. = ^ 2 ) + constant] (positive in 

same sense as a and also relative to thrust azis of the 
rigid aiiT>lane), radians 

rate of change of tail incidence with wing angle of attack 
at constant A 2 (due to fuselage bending under the 
aerodynamic load, 1^.) 

rate of change of tail incidence with wing angle of attack 
at constant 'L^ (due to fuselage bending vinder the loads 
inqposed by the reaction of aft fuselage and empennage 
masses to normal acceleration) 

rate of change of tail incidence with tail angle of attack 
at constant A 2 (due to fuselage bending under the aero- 
dynamic load, L^) 



structural influence coefficient expressing change in tail 
incidence per unit normal acceleration, radians per g 

structural influence coefficient expressing change in tail 
incidence per unit tail load, radians per povmd 


aerodynamic load on horizontal tail (positive when directed 
upward), pounds 

tail length (from airplane center of gravity to aerodynamic 
center of horizontal tail), feet 


free— stream dynamic pressure, pounds per square foot 
tail dynamic pressure, pounds per square foot 
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wing area (including portion covered by fuselage), sqmre feet 
horizontal— tail area, sqmre feet 
true airspeed, feet per second 



tail volume 


cS / 


airplane weight, pomds 

distance from wing aerodynamic center to airplane center of 
gravity (positive when measured forward of center of gravity), 
feet 

spanwise coordinate perpendicxxlar to plane of symmetry, feet 
angle of attack of wing root chord, radians 
angle of attack of tail root chord, radians 
elevator— effectiveness parameter ( | 

lift effectiveness of elevator, per degree 


pitching effectiveness of elevator, per degree 


elevator angle (positive for down deflection), degrees 
downwash angle at the tail, radians 

rate of change of downwash angle at the tail with angle of 
attack 

pitching velocity, radians per second 

nondimens ional spanwise coordinate ( 

\ hl2 

spanwise shift in aerodynamic center position for each wing 

) 

sweep angle of wing quarter-chord line (positive for sweep- 
back), degrees 
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Subscripts 


R rigid airplane 

F flexible airplane 


FUNDAMEMTAL CONSIIEEATIONS 


The purpose of this section of the report is to present the method 
of analysis used in obtaining the calculated results for the example air- 
plane. The material is presented in three main subsections. In the first 
subsection, titled "Stability Equation," the familiar stability equation 
expressing the contribution of wing and tail to static longitudinal sta- 
bility is presented in a modified form to inclxide factors which accovmt 
for the airplane flexibility. In the second subsection, titled "Effect 
of Flexibility on the Stability Parameters," methods are indicated and 
references are given to aid in evaluating the effects of flexibility on 
the factors appearing in the stability equation. In the third subsection 
the effect of flexibility on longitudinal control is discussed. 


Stability Equation 


As is usually done, the present report considers the index of static 
longitxidinal stability to be the partial derivative ^Cj^/dC^. In free 
flight, changes in lift coefficient at constant forward speed must always 
involve a c\irved fli^t path and hence must always involve pitching 
velocity. Further, in the case of a flexible airplane, the normal accel- 
eration associated with a curved flight path will introduce deformations 
due to the loads imposed by the reaction of point masses to normal accel- 
eration. It should be noted that, in general, the effects of such defor- 
mations will be in opposition to the effects associated with deformations 
caused by increase in the aerodynamic loads only. In particular, the 
over-all effect of the mass reactions referred to is to produce loads 
acting normal to the airplane plan form which are distributed over the 
wing and tail spans and along the fuselage. In the discussion to follow^ 
for brevity, these mass effects are referred to as inertial effects but 
should not be confused with the effect of body inertia on dynamic longi- 
tudinal stability where the inertia of the airplane as a whole is con- 
sidered. The sign convention for the present analysis is shown in fig- 
ure 1. 

For a rigid airplane, the stability equati on ■ expressing the contri- 
bution of wing and tail to 30^/30^ is usually written as 
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The additional effect of pitching velocity on Cjj^ results in 


( 1 ) 



V V ^ 


( 2 ) 


If equations (l) and (2) are modified to include effects of wing, tail, 
and fuselage flexihility and are then combined, the following equation 
can he written (see appendix A) 



As can he seen from these equations, flexihility of the airplane struc- 
ture is assumed to affect the following parameters: 

1. Wing aerodynamic center position 

2 . 

3 . 

4 . 

Also, additional parameters are introduced because of fuselage flexibil- 
ity. Of the parameters listed, items 1, 2, and 4 are influenced by wing 
flexibility with item 3 being influenced by tail flexibility. The effect 
of flexibility on these parameters is discussed in detail later in this 
report. 


Wing lift— curve slope 
Tail lift— curve slope 


Rate of change of downwash at the tail 


da 
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It is of some interest to note that the third term of the right- 
hand side of eqmtion (3) and also the additional parameter (3i+/da) 

^t 

in the second term of that equation do not appear when evaluating sta- 
bility for a model in a wind tunnel where the model is physically 
restrained by the support system from developing any accelerations or 
angular velocities. It should be noted, also, that the inertia of vjing 
and tail has not been neglected in developing equation (3), although only 
the effect of fuselage inertia is shown explicitly. Actually, the direct 
effects of wing and tail inertia are considered to be contained implic— 
i'tly in (Cj^ ) and Details of this modification also will 


be referred to later in this report. 


Effect of Flexibility on the Stability Parameters 


Wing aerodynamic center^ .- On a swept-back wing of even moderate 

sweep, wing bending exerts the predominant influence on the aerodynamic 
span load distribution, resulting in an inboard shift in center of load 
for each wing panel. Due to this inboard shift in center of load, the 
aerodynamic-center position for a flexible wing will be ahead of that for 
the rigid wing. As is well known, any forward shift in aerodynamic cen- 
ter will cause a decrease in longitudinal stability. The geometric rela- 
tion between any given spanwise shift and the associated chordwise shift 
along the mean aerodynamic chord is given by the following relation: 


A/^^V-^A^^tanA (4) 

\C J 2 c c/4 ^ 

It is apparent from this expression that the aerodynamic-center shift 
along the chord is most severe at high aspect ratios and high sweep 
angles for a given spanwise shift. A method for determining the span 
load distribution of a flexible swept— back wing, from which the 
aerodynamic-center position of the flexible wing can be determined, is 
given in reference 6. 

Wing lift— curve slope.— On a swept— back wing, the lift— curve slope 
for the flexible wing is usually less than that for a rigid wing. As a 
result, the angle of attack required to reach a given lift coefficient 

^The shift in aerodynamic-center position for the tail is neglected _ 
since it only affects the value of used in the calculation of V. 


8 


MCA. EM A51C19 


is higher for a flexible wing than for a rigid wing. For this reason, 
the tail also is subjected to a hi^er angle of attack (neglecting down- 
wash considerations) so that the stability contribution of a rigid tail 
will be higher on an airplane with a flexible wing than on an airplane 
with a rigid wing. Consequently, the effect of reduction of wing lift- 
curve slope on the stability contribution of the tail is such as to 
increase the stability for an airplane with swept-back wings. 

A method for determining the lift— curve slope for a flexible swept- 
back wing also is given in reference 6. 

Rate of change of downwash at the tail (8c /da).— The rate of change 
of downwash at the tail is dependent upon the span load distribution 
associated with changing angle of attack so that changes in span load 
distribution due to flexibility may have an influence on the average 
downwash at the tail. The rate of chajige of downwash at the tail also 
depends upon the lift— curve slope of the wing, however, so that the over- 
all effect of wing flexibility on cannot be stated even qualita- 

tively for a general case without further analysis. 


Methods for predicting the downwash in the plane of the vortex 
sheet for low lift coefficients for any arbitrary continuous span load 
distribution are given in references 7 and 8. Althou^ the value of 
b^/ba, may veiry considerably in a direction perpendicular to the vortex 
sheet, it was felt that application of such methods to an estimation of 
changes in downwash was at least approximately correct. 

In the method of reference 7 (which is believed to be the simpler), 
the downwash variation be/bCj^ is calculated from a relation which is 
essentially 


8e 

8Cl 


13=1 


(5) 


where the terms a are influence coefficients given in reference 8, 

vn 

and (c,c/CtC„„) are loading coefficients corresponding to given span- 
' 2 ' L av'n 

wise stations as obtained from the span load distribution for the flex- 
ible wing. The span load distribution for the flexible wing can be fo\and 
from the method of reference 6. 


To determine the stability parameter b^/bou, equation (5) is mod- 
ified by introducing (Cj^)^. The res\ilting expression is 
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Tall lift-curve slope.- The effect of flexibility on the lift-c\irve 

slope of a swept tail shows the same qualitative effects as for a swept 
wing. On a swept— back tail, the effect of aeroelasticity will be to 
cause a reduction in lift-curve slope so that the tail lift (and, conse- 
quently, the pitching moment) will be reduced at a given angle of attack. 
The effect of reduction in tail lift-curve slope on the stability of the 
airplane, therefore, is to reduce the stability, whereas the effect of 
reduction in wing lift-curve slope was shown to increase the stability. 

Tail lift— curve slope can be determined by the same method as for 
the wing as given in reference 6. 

Change in tail incidence due to fuselage bending.— From equa- 
tion (3) it can be seen that fuselage flexibility affects the stability 
in general by the introduction of three additional terms giving the 
change in tail incidence due to fuselage bending. These terms are 
(ait/aa)^^, Oi^/aa)^, and (ai^/aa^)^^. Expressions for these terms 

are developed in appendix A but are shown below for convenience. 



(8) 
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Effect of Flexibility on Langitudiml Control 


In the preceding sections of this report, the change in static lon- 
gitudinal stability due to aeroelasticity has been discussed in some 
detail. It is also of interest, however, to consider the effect of flex- 
ibility on elevator angle required for balance = 0). There are two 

effects which must be considered: 

1. The primary effect of the stability changes just discussed. 

2. The secondary effect of elevator deflection in the absence of 

any stability changes. 

The second effect is introduced by the deflection of horizontal tail and 
fuselage under the load produced by elevator deflection. 

The first effect can be evalmted graphically by plotting the sta- 
bility curve (Cjn vs Cl) for the rigid airplane and also the family 
of stability curves for the flexible airplane as calculated at several 
values of dynamic pressure. In turning flight, the elevator angle for 
balance, then, can be calculated from the individual stability curves 
directly. In straight flight, where the dynamic pressure changes with 
lift coefficient, the elevator angle variation can be found from the var- 
iation of Cjjj with Cl given by plotting the flight variation of Cj^ 
with q across the family of curves for the flexible airplane. 

The second effect can be evaluated by means of the following equa- 
tion^ 



^Derivation of this relation is given in appendix B. 


MCA EM A51C19 


11 


In this equation 6g represents the elevator angle required (with flex- 
ibility of the horizontal tail and fuselage present) to maintain the 
same lift on the tail as that produced by an arbitrary elevator deflec- 
tion (&eo^ rigid airplane. The ratio of elevator effectiveness 

for the rigid tail to the elevator effectiveness for the 

flexible tail (SCL^/SSg)^, can be fovind by the relation® 



where k and ACjJ /Ct are constants with magnitude dependent on 

ti tj, 

the structiaral rigidity of the tail. Evaluation of k and ACt* /Ct 

Hi Hp 

(defined in appendix C) involves a knowledge of the span load distribu- 
tion due to elevator deflection and the span load distribution due to 
symmetric twist distribution. Span load distributions for wings and 
flaps of arbitrary plan fonn are given in reference 7 for symmetric flao 
deflection. Span load distributions due to symmetric, continuous twist 
distribution are given in reference 8. The calciaation procedure is 
similar to that contained in reference 6. 

The procedure used to calculate the elevator angle required for 
balance may be simmarized as follows: 

(1) Determine the stability curve for the rigid airplane for 

a given center— of— gravity, position. 

(2) Usi^ the stability given by (l) and the stability change 

given by equation (3), determine the slope of the stabil- 
ity curve for the flexible airplane. Determine the change 
integrating the aeroelastic loadings caused 

by inertia, built-in twist, and camber as found by the 
method of reference 6. With the slope and intercept 
thus determined, plot stability curves for the flexible 
airplane at several values of dynamic pressure. 

3 ■ — - 

Derivation of this relation is given in appendix C. 


12 


MCA EM A51C19 



(3) Plot the flight ■variation of Cj^ •with q across the family 
of curves so oh-tained. 

(^) Determine the ele"vator angle required for balance 5. corre— 
spending to the "variation in C„ with Cy so determined. 

IQ ±j 

(5) Oh'tain the final ele'vator angle required for balance by 
means of equation (lO). 


APPLICATION TO EXAMPLE AIEPIANE 


The method of analysis indicated in the preceding section of the 
report has been applied to an example swep'b-wing airplane known to have 
a relatively flexible structiire. Compressibility considerations in 
regard to the effect on span load distribution were neglected in evaluat- 
ing the aeroelastic effects since a preliminary estimate showed them to 
be of second order (compared to the primary effect of dynamic pressure) 
for the particular configuration studied. The ratio of •tail dynamic 
pressure to free-streeim dynamic pressure was assumed to be equal to 1.0. 
The airplane configuration is shown in figure 2 with the pertinent geo- 
metric parameters indicated. The sweep angles of wing and tail are 35° 
and 33°> respectively; the wing aspect ratio is 9.^3> wing taper ratio 
0.1*2; -tail aspect ratio 4.06, "tail "taper ratio 0.1*23; and "tail volume 
0.672. The effect of the engine nacelles on the aerodynamic span load 
distribution was neglected as was the effect of the fuselage. The elas- 
tic axis for the wing is located at the 38— percen^t— chord line a-nd for 
the tail is located at the 50— percen"t-chord line. The -variation of per- 
tinent s-truct-ural characteristics across the semispan of wing and hori— 
zon-tal tail is presented in figure 3. The structTiral influence coeffi- 
cients associated with fuselage flexibility are defined by the following 


= -0.000031*2, deg/lb 


= 0.1*5, deg/g 
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Effect of Flexibility on Stability 


In the discussion to follow, the material is presented in the fol- 
lowing order: first, an emlimtion of the effects of flexibility on 

each of the parameters appearing in the static— longitudinal-stability 
equation, and second, an evaluation of the over-all effect of flexibil- 
ity for the airplane as a whole. The individual effects are presented 
in figures 4 through 9 aJid are stumuarized in figure 10. The over-all 
effect is presented in figure 11. 

The effect of changes in the stability parameters on the stability 
of the example airplane has been evaluated by ass uming that only the 
parameter under consideration is affected by flexibility. In this way, 
the effect of changes in each parameter can be assessed individually. 

In the discussion of each parameter which follows, effects due to the 
action of aerodynamic loads only are considered first, followed by con- 
sideration of the modifying influence of inertia. 

Wing aerodynamic center.— The shift in wing aerodynamic— center 

position due to wing flexibility is shown in figure 4 together with the 
associated stability change. Inasmuch as the discussion is being 
restricted for the moment to effects due to aerodynamic loads only, 
merely the curves shown for the weightless wing need be considered, since 
these curves are for zero inertial effect. As can be seen, the aerody- 
namic center with inertia absent moves forward from the rigid-^ing value 
of 25 percent of the mean aerodynamic chord until at a dynamic pressure 
of about 650 pomds per sqmre foot the aerodynamic center is at the 
leading edge of the mean aerodynamic chord. At hi^er values of dynamic 
pressTire the aerodynamic center moves even farther ahead. As can be 
seen from the values of stability change, the effect of aerodynamic- 
center shift in itself is very large. For example, at a d/namic pressure 
of 500 pounds per square foot, the neutral point of the wing has shifted 
forward by 20— percent chord, which of itself would introduce a serious 
stability problem. 

It is of some interest to' know how much of the stability change 
with inertia absent is due to bending deformations and how much is due 
to torsional deformations. In order to aid in this comparison, a curve 
has been included in figure 4 showing the contribution of torsion alone. 
As can be seen, torsional deflections are stabilizing, while bending 
deflections are destabilizing. The contribution due to torsion is seen 
to be much smaller than that due to bending. The relative importaince of 
torsion and bending, however, depends on the ratio of torsional to bend- 
ing rigidities and location of the elastic axis, and hence would not 
necessarily be the same for all airplanes. An equally important factor 
to consider is the effect of sweep angle. The extremes of zero sweep 
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and 90° sweep best illustrate the point, since for zero sweep only tor- 
sion is a factor, while for 90° sweep only bending is a factor. 

The effect of inertia on the location of the aerodynamic center 
for the example wing also is included in figure 4 for airplane wing 
loadings of 70 pounds per square foot and 100 pounds per square foot. 

As can be seen from the figure, the effect of wing inertia is only 
mildly alleviating. Although the relieving effect in this case is shown 
to be rather small, the effect for other airplanes may not be of similar 
magnitude,, since the inertia effect depends upon the ratio of wing 
weight to total airplane wei^t in addition to the spanwise weight dis- 
tribution previously mentioned. By reference once more to the case of a 
tailless airplane, it would appear that wing inertia would have a much 
greater relieving effect in that case , since more of the total aiirplane 
weight is in the wings than for conventional airplanes. 


Wing lift-curve slope.- The ratio of flexible to rigid wing lift- 

curve slope and the associated increase in tail contribution to longitu- 
dinal stability is presented in figure 5 as a function of dynamic pres- 
sure. As before, the curves for the weightless wing represent the case 
of zero- inertia effect. At a dynamic pressure of 500 pounds per square 
foot the lift— curve slope with inertia absent is reduced to 64 percent 
of the rigid— wing value. The associated increase in tail stability con- 
tribution amounts to a rearward neutral-point shift of 25 percent. At 
this same dynamic pressure, the stability contribution of the wing 
aerodynamic-center shift (with inertia absent) was shown to be a forward 
shift of 20 percent, or almost the same magnitude, so that the two wing 
factors so far discussed would appear to be largely canceling. Whether 
canceling of these effects will exist in general for all configurations 
cannot be determined at this time. Calcxilations for a fighter configu- 
ration of markedly different geometric and structural characteristics, 
however, resulted in essentially the same relation between these wing 
factors. An interesting extreme to consider is the case of the tailless 
airplane for wMch the second term of the static-longitudinal-stability 
equation does not exist. In this case no canceling of these effects will 
be possible so that the net stability change will be due solely to wing 
aerodynamic— center shift. 

Since the effect on stability of reduction in wing lift— curve slope 
is large with inertia absent, it is of interest to consider the relative 
contribution of bending and torsion. For this purpose a curve is included 
in figure 5 showing the contribution of torsion alone. As can be seen, 
the contribution of torsion to the lift— curve slope is an increase, 
while the larger effect of bending is a decrease. The associated stabil- 
ity changes are shown to be a decrease due to torsion and a much larger 
increase due to bending. 
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The effect of inertia on the wing lift— curve slope also is included 
in figure 5 for airplane wing loadings of 70 and 100 pounds per square 
foot. As can be seen, the effect of wing inertia is much the same as 
was stated for aerodynamic center (i.e., mildly alleviating). 

Bate of change of downwash at the tail.- The variation along the 
swept— tail span of the rate of change of downwash in the plane of the 
vortex sheet with inertial effects absent is presented in figure 6 for 
several values of dynamic pressure. As indicated earlier in this report, 
the curves are based on a method which is applicable only at low lift 
coefficients. The location of the tip of the horizontal tail is indi — 
cated in the figure. As can be seen, large changes in downwash are 
indicated behind the outer sections of the wing and in the plane of 
symmetry; however, the average downwash over the tail is changed only 
slightly. The change in average downwash depends on the ratio of tail 

span to wing span. The downwash factor ^1 - based on the average 

downwash over the tail is presented in figure 7 as a function of dynamic 
pressure along with the associated change in stability contribution of 
the tail. As can be seen from the figure, the change in downwash factor 
is very slight, being of the order of 5 percent at the highest dynamic 
pressure considered. The stability change, as would be expected, is 
correspondingly small and relatively uninrportant compared to the other 
stability factors so far discussed. 

The effect of inertia on downwash at the tail is not shown, since 
the relatively larger downwash changes associated with the aeroelastic 
effects due to aerodynamic loads only were shown to be unimportant. 

Tail lift-curve slope,- The ratio of flexible to rigid tail lift- 

curve slope and the associated decrease in tail stability contribution 
is shown in figure 8 with similar curves for the wing shown for compar- 
ison. As can be seen from the figure, the effect of flexibility on tail 
lift— curve slope is not as pronounced as that for the wing, and, as a 
consequence, the effect on the stability contribution of the tail is 
also correspondingly less. As can also be seen from the figure, the 
inertial effect on the tail is similar to that for the wing. ^ 

5 — 

The figure also indicates the extent to which the inertial effect 
varies with wing loading, since the curve for the weightless wing 
could just as well be labeled w/s = «», A physical explanation here 
is that an airplane of infinite mass will experience zero under 

lifting loads, so that the effect of flexibility for W/s = ® is the 
same as if the airplane were physically restrained and merely pivoted 
about the center of gravity. An airplane having finite mass, however, 
will experience a normal acceleration so that the wing inertial loads 
will affect the wing deflection and hence 
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Tail Incidence change due to fuselage bending.— The 'variation with 

dynamic pressure of an over-all s'tahility term expressing the change in 
tail incidence due to fuselage bending is shown in figure 9- The cur-ves 
showing the contribution of aerodynamic loads only are indicated by 
(Si^/^a)^ = 0. Curves are presented considering the effect of fuselage 

flexibility alone and also including the effects of wing and tail flex- 
ibility. For comparison, ciirves of average downwash and S'tability change 
due to downwash change are also presented. At the hi^er values of 
dynamic pressure the fuselage factor becomes of the same order of magni- 
tude as the average rate of change of downwash at the "tail and therefore 
is seen to be of comparable importance. The effect of including wing 
and tail flexibility in the fuselage factor is to lower the factor 
slightly as shown. By reference to the s'tability c'ur'ves, it can be seen 
that the stability change due to fuselage bending is of much greater 
importance than that due to downwash change, as would be expected from 
the comparison shown on the upper part of the figure. It can also be 
seen that the effect of wing and tail flexibility is to alle-vlate the 
S'tability decrease due to fuselage flexibility. 


As can be seen from the figure, the effect of inertia on the fuse- 
lage factor is very large and consequently of considerable importance. 

It will be remembered that the effect of inertia on the wing and 'tail 
factors was only slight by comparison. It is interesting to note that 
consideration of inertia and of all the flexibilities involved results 
(for the example airplane) in a fuselage factor equal essentially to 
zero even though the aerodynamic contribution is large. In these esti- 
mates of inertial effects, the influence on fuselage factor of wing and 
'tail inertia has been neglected, since these effects are of higher order 
for this airplane. 

Summary of effects.— The effects of wing, tail, and fuselage flex- 
ibility on the longittidinal— S'tability index of the example air- 

plane are summarized in figure 10, showing the impor-tant indi-vidual 
effects which have been discussed. The upper set of cturves presents the 
effects due to aerodynamic load only while the lower set of curves also 
includes the effects due to inertial loads for an airplane wing loading 
of 70 pounds per square foot. As can be seen from the figure, all the 
effects are destabilizing except the effect of reduction in wing lif't— 
curve slope on the s'tability contribution of the -tail. The stability 
changes due to wing aerodynamic— center shift and reduction in wing lift- 
curve slope are shown to be by far the largest effects of those shown. 
Both these results are shown to be true whether inertial effects are 
included or not. 


The over-all effect of flexibility on the static— s'tability index 
SCm/^Ci for the example airplane is shown in figure 11. The curves 
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presented show the changes due to aerodynamic loads only arfl also include 
the effects of inertial loads for an airplane wing loading of 70 pounds 
per square foot. Due to the natvire of the second term of the stability 
equation, the effects shown in figure 10 are not all additive algebrai- 
cally; therefore, the curves of figure 11 were obtained by allowing all 
the factors in the equation to vary simultaneously. As can be seen from 
the figure, inertial effects reduce considerably the stability change 
which would otherwise occur for this airplane. For example, the stabil- 
ity change at a dynamic pressxire of 5OO pounds per square foot is reduced, 
due to inertia, from a neutral— point shift of about I7 percent to about 
7 percent. When the large individml effects are remembered, the over- 
all change due to flexibility for the airplane (as shown by the inertial 
curve) is seen to be relatively small. 

It is of interest to consider the effect of flexibility on the 
fli^t test index (which is associated with ^Ciq/^Cl), since^ 

/^Cm\ 

\ ^89 Jp 

Consideration of the above relation alone leads to the possibility of 
counterbalancing between the effects of aeroelasticity on 
and on therefore, figure 12 is presented to show the over— all 

change in as a function of q for an assumed static margin 

for the rigid airplane of —0.10. Curves are presented for 
SCjn/^be = (^Cm/Sbe)jj and for ^in/Sbe = (^Cin/^6e)j, to show the 
influence of that parameter on A ( ^b^/^Cj^) . As can be seen from the 
figure, the variation of A (3be/^Ci,) with q for the example airplane 
is goveraed primarily by the stability change rather than by loss in 


SThe parameter (^Cj^/^C^) is given by equation (3) of this report. 

The parameter (^Cjj^/^bg) can be found from equation (lO) since, from 

F 

the nature of the analysis leading to that equation, it is clear that 
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control effectiveness. That this conclusion will he true, in genereil, 
for airplanes which suffer loss in stability with increase in q 
perhaps is not immediately obvious, but cem. be seen quite readily from 
recognition of the fact that the intersection of the two curves shown 
in figure 12 at a dynamic pressure of about 600 pounds per square foot 
is due to the flexible airplsune having become neutreilly stable. It is 
apparent that the parameter can have no effect on 

when an airplane has neutral stability since =0. 


Effect of Flexibility on Longitudinal Control 


The effect of flexibility on the elevator single required for 
balance in Ig flight has been evaluated for the example airpleine for a 
static margin of the rigid airplane of O.IO. The semigraphiceil sinalysis 
described earlier was used. The case of ttcmlng flight is not consid- 
ered here since the results for that case are believed to be suffi- 
ciently summarized in figure 12 just discussed. The variation of 
elevator angle with dynamic pressure for steady level flight is pre- 
sented in figure 13(a) for the case where up-elevator only is required 
and in figure 13(b) for the case where both up— an<l downr-elevator are 
required.® Curves are shown for both the rigid and flexible airplane. 

As can be seen from the figure, the effect of aeroelasticity for the 
case of up— elevator only is to increase the up-^levator required over 
the entire speed range by a constant amount of about 0.5°, In the 
case where both up— and downr-elevator are required, the effect of aero- 
elasticity is such as to increase the amount of up— elevator required 
below a d3aiamic pressure of about 400 pounds per square foot ftnti to 
increase the amount of down— elevator required at larger dynamic pres- 
sures. The increase in elevator angle required is less thftTi o.5° 
except at the highest dynamic pressures considered. 

The effect of flexibility of horizontal tail and fuseleige in 
modifying the elevator angle required for balance is shown in figure l4 
as calculated for the example airplane. The effect shown in this figure 
is incliided in the curves of figure 13. The elevatoj>-angle ratio 
increases almost linearly reaching a value of about 1,9 at a dynamic 
pressure of 800 pounds per square foot. 

In all of these calculations, torsion of the elevator has been 
neglected; that is, the surface is assumed to be infinitely rigid to 
hinge moments. 

®The determining factor in each case is the value of for the 

rigid airplane. In case 1, was assumed eq^ial to zero, while 

for case 2, CnQ was assumed eqrial to +0.03. 
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CONCLUDING EEMAEEB 


The effect of flexibility on static longitudinal stability and 
control of an airplane is such as to preclude einy real generalizations 
of the results presented herein. Although the over-all aeroelastic 
effect on stability for the example airplane was foiind to be smll coia- 
pared to the individual effects, it cannot be said that like calcula- 
tions for any airplane will also yield a small effect. It can be said, 
however, that, for any practical swept-wing airplane with a tail, the 
stability change due to a shift in wing aerodynamic center will be 
destebilizing, while the change in wing lift-curve slope will be stabi- 
lizing, so that a certain amount of counterbalancing between these major 
effects will always be present. As can be seen from the simple stabil— 
ity equation en 5 >loyed in this analysis, the degree of completeness of 
the coTinterbalancing depends directly on the size, plan form, and loca- 
tion of the tail as they affect the stability contributed by the tail. 
Therefore, in the design of airplanes for which wing flexibility would 
be expected to exert a large influence on stability, it would appear 
that a minimum over-all aeroelastic effect may be accomplished more 
advantageously by design changes to the horizontal tail than by such 
changes to the wing. 

The calcinations presented herein with regard to the effect of flex- 
ibility on longitudinal control show the effect to be of little real 
in 5 )ortance for the exaniple airplane in Ig flight, although in turning 
flight the effect is shown to be coiamensurate with the stability loss, 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Moffett Field, Calif. 
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APPENDIX A 

EYALUATION OF PAEAMETERS ASSOCIATED WITH 
FUSELAGE FLEXIBILITT 


In order to show how the parameters associated with fuselage flex— 
ihility enter the longitudinal— stahility equation, equation (3) of the 
text will he derived here. If the pitching-moment contribution of wing 
and tail are considered in nondimens icaaal form (i.e., in terms of Cj^, 
Cl^, etc.), the pitching-moment coefficient of the comhimtion can he 
written as : 


Since a flexible airplane is being considered, (C^) and (Cj^) 

F F 

should be expressed in terms of parameters applying to a flexible struc- 
ture, so that 

■ 

('■'i ■ 

where 

CL^ = 0^ - 


^ ^F 

(A2) 

F 

(A3) 

^F + it + ^t 1 ' 

(a4) 


In eqmtion (a 4), it should be recognized that the downwash angle Cj. 
may be affected by the changes in span load distribution associated with 
flexibility and also that the tail incidence i^ is affected by fuselage 
flexibility. Substituting eqmtion (a 4) into eqmtion (A3) followed by 
substitution of the modified eqmtion (A3) into eqmtion (Al), 
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wiliti r©sp6cb to f*8.ctoriiig ^ f^rom tii6 

(X TJT 


second term, equation (A5) becomes 

.) 

:/ 5; 


5£m A _ V H 

\ ^Cl y V c 4 ( c-r. ^ 


'Y 


1 - ^ + ^ + ii M. 

V^a/ Sap V Sop J 

F 


v3i (a6) 


Before proceeding further, it is necessary to evaluate Si^/Sa. To do 
this, the tail incidence (i-t^ written 


^ \ SLt / 




(A7) 


where (Si^/SLj.) and (Si+/SA 2 ;)^ are structural influence coeffi- 
A 2 

cients associated with fuselage bending. Differentiating equation (AT) 
with respect to Op, 


dAr 


dop \ SL(. / dop \ ^A^ / dop 

^t 


(a8) 


It should he noted here that, although equation (a 6) contains the par- 
tial derivative Si^/Sop, the total derivative of equation (AT) must he 

taken because of the dependence of e , i . , and 0 on a in the 
expression for L^. which can he written as 


Gv ) + it + Y 


(A9) 


V 

The derivative of equation (A9) with respect to a^, is then 


— = (cr ) Fi ) 

da^, \ itc, . / L \ da / 


+ iit + It de’ 


dOp 


dOp J 


‘It^t 


(AlO) 
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For notational consistency with equation (AT), de/dop and dd/dop can 
he written as de/3ap and since at constant q the quanti- 

ties e and § are functions of Op only. With these substitutions 
in equation (AlO) and knowing that 



(All) 


from differentiation of 




(A12) 


equation (a 8) can be written after some rearreuigement as 



For notational convenience in expressing t^ final equation, let 
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and 





so that eqijation (A14) can he written as 


+(^') 

dop \ 3a / \ 3o+ / V aa \ / 

2i ^ 


(Al4) 


By substitution of equation (ai 4) into equation (a 6) and after rearrange- 
ment, a final expression for aCjjj/aCj^ may be written as 


©. ■« - 




F 


T7* 


V 


<lt 




q acp 


(A15) 
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APPENDIX B 

EVALUATION OF ELEVATOE ANGLE EEQUIRED FOE BALANCE 


In. evalnating the elevator angle required for balance (Cj^ = 0) 
for a flexible airplane, it is necessary to consider the effect that 
elevator deflection will have on aeroelastic distortion of the horizon- 
tal tail and fuselage since such distortion will change the elevator 
angle required to balance the pitching moment existing with elevator 
neutral. If the elevator is assumed to be deflected by an arbitrary 
angle 5^ with the fuselage and horizontal tail fixed in position, the 

lift coefficient on the tail due to elevator deflection may be written 
as 



(Bl) 


If the horizontal tail is now allowed to relax (with the fuselage still 
position), the lift on the tail will change due to distortion 
of the structure. The elevator angle required to maintain the same lift 
on the tail as that given by equation (Bl) is defined by 



If the fuselage now is allowed to relax, the lift on the tail will change 
due to a change in tail incidence. The additional elevator angle required 
to maintain the lift at the value given by equation (Bl) can be written 
as 



(b4) 
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where 



The final elevator angle required is then 


^e ” ^ei ^62 


(b6) 


By suhstitution in equation (b 6), the following expression is obtained 
for the elevator angle required on a flexible airplane to maintain the 
same lift on the tail produced by an arbitrary elevator deflection on a 
rigid airplane. 
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APPENDIX C 

EVAIUATION OF LIFT EFFECTIVENESS FOR FLEXIBLE TAIL 


The lift effectiveness of an elevator (3Cj^/^5g) for the case of 

a flexible horizontal tail will differ from that for a rigid tail due to 
distortion of the struct^e.® The streamwise twist of the structure can 
be found as in reference 6 by applying a relaxation approach to the prob- 
lem. Using this approach, the twist distribution for a swept-back tail 
can be written in series form as 


€(Tl) = +A 62 (ti) + . . . (Cl) 

where 

e(t^) twist of the flexible wing 

€o(t)) twist produced by the loading for the rigid wing 
A€i('h) twist produced by the loading associated with 
Ae 2 (Tj) twist produced by the loading associated with 
etc. 

Since g(ti) can be written in series form, it is apparent that a simi- 
lar expression can be written for the lift coefficient produced by a 
given elevator deflection for the flexible wing, so that 


Ct = Ct - A[)t + ACt . . . 

\ \ ■ ti H 2 


(C2) 


where 

C-r lift coefficient for the flexible tail 

S’ 

Ct lift coefficient for the rigid tail 

\ 

AC lift coefficient obtained by integrating the loading associated 

with g^(t]) 

^Elevator distortion is neglected in this analysis. 
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^L+ lift coefficient obtained by integrating the loading associated 
^ with Aei(Ti) 

etc. 

If the terms of equation (Cl) are related by a constant of proportion- 
ality (which is the usual case),"^ equation (C2) can be written as 



since the series 1 — K + K ... is merely an expansion of l/l+K. 
The ratio of the lift effectiveness for the flexible tail to that for 
the rigid tail can be obtained from equation (C^^) by dividing by Cl 
T he equation becomes ^ 






1 

1 +K 


(C5) 


7 

The line of load application for the first term is the centroid 

of the chordwise loading produced by elevator deflection while the 
line of load application for the remaining terms is the aerodynamic 
center of the section. The relative contribution of bending and tor- 
sion, therefore, is not the same in €q(tj) as in succeeding terms. 
Calculations made to date, however, have shown the constant of pro- 
portionality to apply to the relation between ^q(ii) and Aei(q) 

as well as to the terms beyond Aej^(q). There may be some configura- 
tions for which the proportionality will be limited to Aei(q) and 
succeeding terms, however. 
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Iqmtion (C5) can "be simplified since K and 
to and ‘becomes 



are proportional 



where 

tail. 



and k are evaluated for unit dynamic pressure at the 


MCA BM A51C19 


29 


N 


HEFEREIOIS 


1. Collar, A. E. : The Expanding Domain of Aeroelasticity. British 

A.R.C. 9701, S.& C. 2022. June 11, 19k6. 

2. Collar, A. E., and Irinsted, F. : The Effects of Structural Flex- 

ibility of Tailplane, Elevator, and Fuselage on longitudinal 
Control and Stability. E.&, M. 2010, British A. E.C. , September 
191^2. 

3. Lyon, H. M. : A Method of Estimating the Effect of Aeroeiastic 

Distortion of a Sweptback Wing on Stability and Control Deriva- 
tives. British E. A. E. T.N. Aero. l8l3, July 1946. 

4. Lyon, H.M. , and Eipley, J. : A General Survey of the Effects of 

Flexibility of the Fuselage, Tail Unit, and Control Systems on 
Longitudinal Stability and Control. Eep. No. Aero. 2065, R.A.E. 
British, July 1945. 

5. Pugsley, A. G. : Influence of Wing Elasticity Upon Longitudinal 

Stability of an Aeroplane. E.&.M. 1548, British A. E.C. , 1933. 

6. Skoog, R. B., and Brown, H. H. : A Method for the Determination 

of the Span Load Distribution of a Flexible Swept Wing at Subsonic 
Speeds. NAGA TN 2222, I 95 I. 

7. DeYoung, John: Theoretical Symmetric Span Loading Due to Flap 

Deflection for Wings of Arbitrary Plan Form at Subsonic Speeds. 
HACA TN 2278, 1951 . 

8; Diederich, F. W. : Charts and Tables for Use in Calculations of 

Dovmwash of Winigs of Arbitrecry Plan Form, NA.CA TN 2353^ 1951. 

9* DeYoung, John, and Harper, C, W, : TheoreticeCL Symmetric Span 

Loadiig at Subsonic Speeds for Wings Having Arbitrary Plan Form, 
NACA Eep, 921, 1948, (Formerly NACA TN»s 1476, 1491, and 1772) 



Page Intentionally left blank 


Page intentionally left blank 






Figure 2.- Geometric characteristics of example airplane. 


Dead -weight shear, tb Torsional rigidity, GJ, Bending rigidity, El, 
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Spanwise station, 
(a) Wing. 



Figure 3.- Variation in pertinent structural character- 
istics across the semispan for the wing and hori- 
zontal tail of the example airplane. 
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Figure 4 -Variation of wing aerodynamic center with dynamic pressure for the example 

airplane . 




Torsion only (weightless wing) 



Figure 5. - Variation of the ratio of flexible to rigid wing lift - curve slope with 

dynamic pressure for the example airplane. ■ 
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Figure 6 - Variation along the swept- tail span of the rate of change of down wash in 
the plane of the vortex sheet for several values of dynamic pressure for the 
example airplane as affected by aerodynamic loads only. 
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Figure 8- Variation of the ratio of flexible to rigid tail lift -curve slope and 
associated stability change with dynamic pressure for the example airplane^ 
with similiar curves for the wina shown for comoarison . 
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Dynamic pressure, q, !b/sq ft 

Figure 9 - Variation with dynamic pressure of the factors associated with 
fuselage flexibility and the stability change due to these factors for 
the example airplane. 
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Figure iO- Summary of the individual effects of the various 
parameters involved on stability of the example airplane 
as affected by aerodynamic toads only and also as af- 
fected by both aerodynamic and inertial toads. 
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pure 12.- Over- all effect of flexibility on the rate of change of elevator angle with lift 
coefficient in turning flight for the example airplane assuming static margin for the 
rigid airplane of -OJO 
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